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ABSTRACT 

The naked-eye star 55 Cancri hosts a planetary system with five known planets, including a hot super-Earth (55 Cnc e) extremely close 
to its star and a farther out giant planet (55 Cnc b), found in milder irradiation conditions with respect to other known hot Jupiters. This 
system raises important questions on the evolution of atmospheres for close-in exoplanets, and the dependence with planetary mass 
and irradiation. These questions can be addressed by Lyman-o' transit observations of the extended hydrogen planetary atmospheres, 
complemented by contemporaneous measurements of the stellar X-ray flux. In fact, planet 'e' has been detected in transit, suggesting 
the system is seen nearly edge-on. Yet, planet 'b' has not been observed in transit so far. Here, we report on Hubble Space Telescope 
STIS Lya and Chandra ACIS-S X-ray observations of 55 Cnc. These simultaneous observations cover two transits of 55 Cnc e and 
two inferior conjunctions of 55 Cnc b. They reveal the star as a bright Lya target and a variable X-ray source. While no significant 
signal is detected during the transits of 55 Cnc e, we detect a surprising Lya absorption of 7.5 ± 1.8% (4.2 a) at inferior conjunctions 
of 55 Cnc b. The absorption is only detected over the range of Doppler velocities where the stellar radiation repels hydrogen atoms 
towards the observer. We calculate a false-alarm probability of 4.4%, which takes into account the a-priori unknown transit parameters. 
This result suggests the possibility that 55 Cnc b has an extended upper H i atmosphere, which undergoes partial transits when the 
planet grazes the stellar disc. If confirmed, it would show that planets cooler than hot Jupiters can also have extended atmospheres. 

Key words, planetary systems - stars: individual: 55 Cancri - planets and satellites: atmospheres - techniques: spectroscopic 



1. Introduction 

A large fraction (~27%) of exoplanets are found at close separa- 
tions(0.01-0. 1 au) from their host stars, as a result of a detection 
biasPjThese hot planets sustain high levels of extreme ultraviolet 
(EUV) and X-ray radiations, which can impact their atmospheric 
properties, dynamics, and ultimately their fate. Observations 
of transits in the stellar Lyman-or (Lya) emission line of neu- 
tral hydrogen (Hi) have revealed extended upper atmospheres 
on the hot giant exoplanets HP 209458b ([Vidal-Madjar et al.| 



2003 , 2004 2008 ; Ben-Jaffel 2007 2008 ; Ehrenreich et al. 2008 ) 
and HD 189733b ( [Lecavelier des Etangs et al.||2010| |2012) . 
These extended H i envelopes are also observed along with heav- 
ier atoms and ions (| Vidal-Madjar et al.||2004} |Ben-Jaffel &| 
|Sona Hosseini|2010 Linsky et al.|2010 1, also possibly detected 
on Wasp-12b (Fossati et al. 2010J. The escape of heavy ele- 
ments suggests the upper atmosp heres are in a hydrodynamica l 
'blowoffPfstate (|Garcia Munoz|2007| [Murray-Clay et al.|2009 1. 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are asso- 
ciated with GO/DD programme #12681. 

1 The velocimetric and photometric detection techniques are more 
efficient for short-period planets. 

2 A discussion about the possible meanings of 'blowoff ' can be found 
in|Ti an et al"] ( |2008) 



Such hydrodynamic escape is thought to have occurred for much 
cooler planets in the early Solar System (Lam mer et al.| [2008). 
This raises two interesting questions: (i) Is there a critical dis- 
tance to the star below which an atmosphere would experience 
a sustained hydrodynamical escape (Koski nen et al.|2 007)? (ii) 
Is there a critical mass threshold below which a planet can be 
totally stripped of its atmosphere (Lecavelier des Etangs et al. 
[20041 |Owen & Jackson|2012) i? 

The 55 Cancri planetary system can help address both ques- 
tions. This star (G8V) is n earby (12.3 pc) , bright (V = 5.95), and 
a host to five exoplanets ( Fischer et al.||2008 l. Planet 'e' is the 
closest to the star, a hot super-Earth (M e = 7.81^S - i| M e ) at only 
0.0156 au (Dawson & Fabrycky 2010), recently detected in tran- 
sit fWinn et al.|2011||Demory et al.|20"T2| l. Planet 'b' is the sec- 
ond closest, a Jupiter-mass object (Mj, = 0.8 M^) at 0.115 au, 
close to the predicted limit of 0. 15 au where 'a sharp breakdown 
in atmospheric stability' should occur for jovian planets around 
solar-type stars (Koskin en et al.|2007| l. Studying the exospheric 
properties of 55 Cnc e and b could therefore provide key insights 
into the critical mass and distance at which atmospheric escape 
may dramatically impact the atmospheric evolution of exoplan- 
ets. However, 55 Cnc b has not been detected in transit so far. 

In this article, we report on Hubble Space Telescope (HST) 
Lya observations of transits of 55 Cnc e and inferior conjunc- 
tions of 55 Cnc b (Sect. They are complemented by con- 
temporaneous Chandra observations aimed at monitoring the X- 
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ray irradiation received by the planets. While the Lya observa- 
tions show no evidence for a transit signal during the transits of 
planet e, they hint at a transit-like signal detected in phase with 
inferior conjunctions of planet b (Sect. [3]), suggesting that the ex- 
tended upper atmosphere of 55 Cnc b undergoes partial transits 
(Sect.@}. 

2. Observations and data analysis 

2.1. Hubble Space Telescope Lya observations 

We obtained HST time (GO/DD 12681) to observe 55 Cnc dur- 
ing two visits on 2012-03-07 and 201 2-04-05 with the Spac e 
Telescope Imaging Spectrograph (STIS; Woodg ate et al.|1998] >. 
Each visit consists of four consecutive HST orbits. The data 
were recorded in time-tag mode with exposure times of 2040 
to 3 005 s on the Far Ultraviolet Multi-Anode MicroChannel 
Array (FUV-MAMA) detector. Incoming light was diffracted 
with the long 52" x 071 slit and the G140M grating, yielding 
first-order spectra centered on 1 222 A and ranging from 1 194 
to 1 249 A. The medium spectral resolution obtained is ~ 10000 
(~ 3010ns- 1 ). 

Figure [T] shows an example of a STIS/G140M spectrum of 
55 Cnc. The Lya line (1 215.67 A) is, by far, the most prominent 
feature in the whole emission spectrum. Between 1 214.5 and 
1 217 A, we measured mean Lya fluxes of (4.19 ± 0.02) x 10~ 12 
and (4.36 + 0.02) x 10~ 12 erg cirT 2 s _1 in visits 1 and 2, re- 
spectively. The quoted uncertainties are provided by the STIS 
pipeline. They are not representative of the observed flux scat- 
ter, as shown in Fig.|2]and discussed in the next paragraph. These 
Lya flux levels are brighter than for any other known transit- 
hosting stars, as shown in Fig. [3] The Lya feature has a typical 
double peak due to absorption in the line centre by the inter- 
stellar H i along the line of sight. Here, the interstellar absorp- 
tion is narrow, and in particular the 'blue peak' of the line is 
much less absorbed by the interstellar medium than in the case 
of HD 189733, allowing us to probe the line at lower Doppler 
velocities. Other stellar emission from ionised silicon (Si m at 
1 206.5 A), oxygen (O v at 1 218.3 A), and nitrogen (N v doublet 
at 1 238.8 and 1 242.8 A) are also detected. These lines do not 
show any significant variations during the observations and in 
the following, we will dedicate our attention to the Lya line. 

Taking advantage of the time-tag data acquisition mode, each 
orbit-long exposure was split into 500-s sub-exposures, for a to- 
tal of 22 sub-exposures per visit. This enabled a monitoring of 
the Lya flux during each HST orbit (Fig. [2}, which revealed a 
systematic effect clearly seen in visit 1: the Lya flux increases 
with HST orbital phase (Fig. [4}, introducing a scatter of 7% and 
3% around the average fluxes quoted above for visits 1 and 2, 
respectively. This effect (of a maximal amplitude of ~ 20%) is 
attributed to optical telescope assembly 'breathing' - secondary 
mirror motions due to orbital thermal variations (Kimbl e et al.| 
|1998| . We determined that this effect can be modelled with a 
quadratic function of the orbital phase (see also Sing et al. 201 1) . 
The fit was performed on a sample trimmed from the exposures 
obtained near the inferior conjunction of 55 Cnc b, so as to 
avoid a possible signal contamination, yielding a^ 2 of 3.85 for 
v = 8 - 3 degrees of freedom (y 2 = 0.77). 

In contrast, visit 2 is not impacted by telescope breathing. We 
checked that fitting linear or quadratic functions to the phase- 
folded Lya light curve (also trimmed from values close to the 
inferior conjunction of planet 'b') does not statistically improve 
the dispersion measured in the 'raw' light curve. Correction fac- 
tors are derived as a function of the HST orbital phase, as the ra- 
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Fig. 1. STIS/G140M spectrum of 55 Cnc. 
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Fig. 2. Top panels: Lyman-a flux integrated over the whole Lya 
line as a function of time for 500-s sub-exposures in visits 1 
(circles, left) and 2 (squares, right). The different colours mark 
the different HST orbits. Bottom panels: Lyman-a flux corrected 
for the 'breathing effect' . Dotted lines show the mean fluxes in 
visits 1 and 2. 



tio between the polynomial value to the mean polynomial value 
over the observed phase range. The detrended Lya light curves 
extracted from the corrected spectra are shown in Fig. [2] 

2.2. Chandra X-ray observations 

In order to monitor the stellar X-ray radiation, we obtained 
Chandra ACIS-S coverage ( Weisskopf et al. 2002 Garmire et al. 
[2003] ) of each HST visit (Chandra ObsIDs 14401 & 14402). The 
observations were scheduled to begin 4 hours before mid-transit 
of 55 Cnc e, with precise start times of 2012-03-07 UT 02:13 
and 2012-04-05 UT 12:50. The scheduled exposure times were 
20 ks, but the first observation was cut short to 12 ks by safing of 
ACIS due to high Solar radiation levels. Since 55 Cnc is a bright 
optical source, we operated the ACIS-S CCD in very faint (VF) 
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The background was estimated from a large adjacent region on 
the same CCD selected to be free from other sources. The count 
rate was found to be 0.0187 + 0.0015 s _1 . 

3. Results 

3.1. 55 Cancrib 
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Fig. 3. Comparison of STIS/Lyo- spectra of stars hosting tran- 
siting exoplanets: the G8V star 55 Cnc at 12.3 pc (this work, 



blue), the K1-K2 star HD 189733 at 19.3 pc ( |Lecavelier des 
Etangs et al."1|20T2] oran ge), the M3.5V star GJ 436 at 10.2 pc 



(Ehrenrei ch et al.||2011[ red), and the GOV star HD 209458 at 
49.6 pc ( [Vidal-Madjar et al.|2003| yellow). 
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Fig. 4. Lyman-ff flux phase-folded on the HST orbital period. 
Filled symbols are exposures taken near the inferior conjunction 
of 55 Cnc b, which are excluded from the fit. The dashed lines 
indicate the best-fit quadratic function for visit 1 ('breathing ef- 
fect'; left) and constant for visit 2 (right). 



mode and selected a 1/8 sub-array in order to minimise optical 
loading. 

We analysed the Chandra data using CIAO version 4.4. A 
weak soft X-ray source was detected at the proper-motion cor- 
rected co-ordinates of 55 Cnc in both observations. We extracted 
source counts using a region of radius I'.'l, and the background 
was estimated using two large circular regions of radius 27" ad- 
jacent to the target and aligned along the sub-array. The source 
had consistent background subtracted count rates in the two ob- 
servations of 0.0019 ± 0.0005 and 0.0022 ± 0.0004 s _1 . Chandra 
ACIS-S light curves of 55 Cnc are presented in Figs. |5]&[8] 

For comparison, we also analysed an archival XMM-Newton 
X-ray observation of 55 Cnc from April 2 009 (XMM-Newton 
ObsID. 0551020801) previously reported by|Sanz-Forcada et al. 



( |201 l| l. We used SAS version 12.0 and extracted source counts 
from the EPIC-pn camera using a circular region of radius 19V 2. 



The orbital properties o f 55 Cnc b (|Butler et al.|1997| l were ac- 
curately constrained by Fisc her et al.| ( |2008[ ), from 18 years of 
velocimetric measurements (from 1989 to 2007) obtained at the 
Lick and Keck observatories. These authors ruled out central 
transits of planets b, c, and d. |Dawson & Fabrycky ( |2010| l re- 
analysed the data and revised the orbital period of 55 Cnc e, from 
2.8 to 0.74 days. Transits of planet e have been subsequently de- 
tected (Winn et al. 201 1| |Demory et al. 2012| l, suggesting the 
whole system is seen nearly edge-on. There is thus a possibility 
that the second closest planet, 55 Cnc b, could graze the stellar 
disc at its closest projected approach - the inferior conjunction. 

3.1.1. Constraints on the inferior conjunction times 



Dawson & Fabrycky (2010 1 provide an ephemeris for 55 Cnc b, 
predicting an inferior conjunction at BJD 2453 092.752 ± 
0.023 (2004-03-28). To update their prediction as accurately 
as possible, we pooled together Keck and Lick ( |Fischer et al.| 
[2008] ), HET HMcArthur et al.|2004| >, and SOPHIE radial veloc- 
ities (Moutou et al., in preparation). We performed a 5-planet 
Keplerian fit to these data and introduced a systematic 'jitter' 
for each data set, which is free to vary and quadratically added 
to the quoted uncertainties on the radial velocities, so to weigh 
the relative contribution of each time series according to its 
true precision. We obtained for 55 Cnc b an orbital period of 
14.651 364 + 0.000074 days, and predict inferior conjunctions 
during HST visits on BJD 2455 993.770 ± 0.013 (2012-03-07 
UT 06:28+19 min) and BJD 2456023.073 ± 0.013 (2012-04-05 
UT 13:44+19 min). 

Combining both HST visits provides an excellent coverage 
of the inferior conjunction of 55 Cnc b. As for HD 209458b and 
HD 189733b, we would expect any transit signatures as absorp- 
tion over a limited and blueshifted range of velocities in the Lya 
line. We scanned the spectra in search for transit-like signatures 
in the vicinity of the inferior conjunction. For each wavelength 
range, we calculated a phase-folded transit light curve. Since 
the average Lya flux level varies between the two epochs (see 
Fig. [2J, the fluxes are normalised in each visit by the mean 'out- 
of-transit' flux (which is determined a posteriori). At the pre- 
dicted time of the inferior conjunction, we found a blueshifted 
absorption signal of 7.5+1 .8% between 1 215.36 and 1 215.67 A, 
or between Doppler velocities of -76.5 and km s _1 , as can be 
seen on the phase-folded light curve (Fig. [5} and on the trans- 
mission spectrum (Fig|6]l. 

3.1 .2. Impact of background subtraction 

This wavelength range is contaminated by the terrestrial air 
glow, which shows large - yet reproducible - variations during 
the HST orbits within a visit. The air glow impact on Lya tran- 
sits has been extensively discussed in the case of HD 209458 
~ 20581 |Ben-Jaffel||2007) . However, unlike 



(Vidal-Madjar et al. 



HD 209458 where the Lya line is about as bright as the air 
glow, the Lya emission of 55 Cnc is ~ 20 times brighter than 
the air glow. We nevertheless performed a battery of tests to as- 
sess the possible impact of the air glow subtraction to the re- 
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duced spectra. The STIS reduction pipeline takes advantage of 
the spatial extent of the air glow on the detector y-axis to sub- 
tract it, together with the background. The background level is 
estimated in a region located +300 pixels away from the cen- 
tre of the spectral extraction region (Bostr oem & Proffitt|201 l| l. 
Ben-Jaffel ( 2007 ) called attention to the possible flux variations 
on the detector along the y-axis, which represent - according to 
this author- an incompressible uncertainty of 5% on the flux 
level of extended sources. We mapped the air glow along the y- 
axis for each HST orbit and found that systematic variations are 
indeed present. To quantify their effect, we ran several data re- 
ductions, extracting the background at different locations along 
the detector y axis. These tests show that the air glow - and more 
generally the background - subtraction (i) does not impact the 
detection of a transit signal for 55 Cnc b, yet (ii) induces a mean 
dispersion of ~ 3% on the flux measured on each 500-s sub- 
exposure. This dispersion was quadratically added to the flux 
uncertainties, and propagated throughout our calculations. 



3.1.3. Transit-like signal and false-alarm probability 

In order to estimate the mean 'in-transit' absorption, we need to 
define the transit contacts, i.e., which exposures are 'in transit' 
and which ones are 'out of transit'. Because the transit is not 
detected in the optical and the Lya transit duration is not known, 
a systematic exploration of possible start and end times of the 
transit has been undertaken. The largest absorption signal (7.5 ± 
1.8%) is obtained from ~ 2.3 h before to ~ 1 h after the inferior 
conjunction of 55 Cnc b. 

False-alarm probabilities (FAPs) should be associated with 
the measured signal, to take the choices of a-priori unknown pa- 
rameters into account. These choices are (i) the transit start (Ti b ) 
and end (Tiy b ) times, and (ii) the wavelength domain in the Lya 
line. The FAP associated with (i) is calculated with a boot-strap 
method, randomising the choices of Tj b and Ti Vb . In each trial, 
the transit depth is estimated and the associated signal-to-noise 
ratio is compared to the fiducial value (S /N = 4.2). After 10 5 
trials, we found a FAP of 1.8%. Next, to estimate the FAP asso- 
ciated with (ii), we generated new spectra for all of the 500-s ex- 
posures, combining a reference averaged spectrum with random 
noise. Then, we searched the blue part of the Lya line in gener- 
ated spectra for a wavelength region where the in-transit absorp- 
tion (that would be due here to noise flucutation) would be more 
significant than the fiducial signal. For 5 x 10 4 trials, we found a 
FAP of 2.6%. Since these FAPs are independent, the final prob- 
ability that the signal is real is (1 - 0.018) x (1 - 0.026) = 95.6% 
(FAP of 4.4%). 



3.2. 55 Cancri e 

A similar analysis is performed for 55 Cnc e. Figure [7] shows 
the absorption depth 6=1- F; n /F out as a function of wave- 
length (the transmission spectrum) for planet e. The 'best' sig- 
nal, 7.0 ± 5.3%, is found between 1 215.56 and 1 215.67 A or 
between -26.6 and km s _1 , with an associated FAP of 89%. 
The weak significance of this signal coupled with the high FAP 
call for a spurious origin. 

As a test, we calculated the transit signal of 55 Cnc e over 
the same wavelength interval used for 55 Cnc b, i.e., between 
1 215.36 and 1 215.67 A. On this interval, where it is possible to 
measure a transit depth for 55 Cnc e with the same precision as 
for 55 Cnc b, we obtain an absorption depth of 0.3 ± 2.4%. The 
corresponding light curve is plotted in Fig.[8](top panel). 
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Fig. 5. Top. Lyman-a light curve integrated between 1 215.36 
and 1215.67 A or between -76.5 and km s _1 . The expo- 
sures from visits 1 (small circles) and 2 (small squares) have 
been phase-folded around the period of 55 Cnc b. The different 
colours (grey, blue, green, red) correspond to the four orbits in 
each visit. Exposures considered in and out of transit are dif- 
ferenciated by filled and empty symbols, respectively. The large 
black symbols (circles for visit 1, squares for visit 2) are the 
flux averages per orbit. The sky-blue hatched region shows the 
±19-min uncertainty on the inferior conjunction time. Bottom. 
Chandra X-ray count rate integrated over ~ 1-h bins. 
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Fig. 7. Absorption depths during transits of 55 Cnc e, as a func- 
tion of wavelength along the Lya line. 



Our conclusion is that we do not detect any transit signal 
from 55 Cnc e at Lya. 



3.3. 55 Cancri X-ray emission 

Inspection of the Chandra ACIS-S spectra of 55 Cnc shows 
it to be a very soft source (as expected for a relatively inac- 
tive late-type star) with the majority of detected events lying 
softer than the limit of the Chandra ACIS-S spectral calibration 
(0.243 keV). Analysing the spectrum from the longer second ob- 
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Fig. 6. Top. Lyman-a line of 55 Cnc. Weighted co-addition of exposures flagged 'in' (F[ n ; solid blue line) and 'out' (F out ; solid 
red line) of the possible transit of 55 Cnc b's extended atmosphere. The solid black lines show the stellar intrinsic emission profile 
before (thin) and after (thick) absorption by the interstellar medium, in terms of flux (left axis) or/3 ratio between radiation pressure 
and gravity (right axis). The absorbed profile is convoluted with the instrumental line spread function to yield the green profile, 
which is fit to the data. Dashed vertical lines delimit the region over which the flux is integrated to produce the light curve in Fig. [5] 
Bottom. Transmission spectrum showing the absorption depth 6=1- Fi n /F out during the possible transit of 55 Cnc b. 



servation (using XSPEC version 12.7), we found a marginally 
acceptable fit using a single-temperature collisionally-ionised 
plasma model (APEC; Smith et al. 2001). Given the low num- 
ber of counts we minimised using the Cash statistic and tested 
goodness of fit using Monte-Carlo simulations. For a best fitting 
single temperature of 0.09 ± 0.02 keV we found 80% of sim- 
ulated spectra yield lower fit statistics. The [0.25, 2.5]-keV or 
[5, 50]-A flux was 2.4 x 10~ 14 erg s" 1 cm" " 
X-ray luminosity in this band of 4 x 10 26 



, corresponding to an 
erg s . 



In order to make the best comparison with the X-ray flux at 
the time of the XMM-Newton observation in April 2009, we car- 
ried out a simultaneous fit of the XMM-Newton EPIC-pn spec- 
trum and our second (longer) Chandra ACIS-S observation. We 
used a two-temperature APEC model, forcing the temperatures 
to be the same for the two observations but allowing the nor- 
malisations vary indepdendently. It was immediately obvious 
that the X-ray emission at the time of the Chandra observation 
was both softer and fainter than at the time of the XMM-Newton 
observation. An acceptable fit was found (with 62% of Monte- 
Carlo realisations of the model spectrum yielding lower values 
of the Cash statistic) with best fitting temperatures of 0.09 + 0.01 



and 0.6 + 0.1 keV. The normalisation of the low-temperature 
component was a factor 1.7 lower for the Chandra data com- 
pared to the XMM-Newton data, while the normalisation of the 
high-temperature component was a factor of 45 lower. This sug- 
gests that the EUV emission was relatively similar at the times 
of the two observations, while the X-ray flux was much lower 
at the time of the Chandra observation. Averaged across the 
[0.25, 2.5]-keV ([5,50]-A) band, the X-ray emission was a fac- 
tor 2.4 lower at the time of our Chandra observations than at the 
time of the XMM-Newton observation in April 2009. 

It is clear that 55 Cnc is a variable EUV and X-ray source, 
with weaker irradiation of 55 Cnc e and b in March/April 2012 
than in April 2009. The X-ray flux was consistent between the 
two HST visits, however, and no flares were detected, giving us 
confidence that it is reasonable to assume that the mass loss rates 
from the planets were also similar at the two HST epochs. 

Together with the detection of variable mass loss from 
HD 189733b ( Lecavelier des Etangs et aT7||2012| >, our results 
underline the need for future Lyman-a studies of extended 
planetary atmospheres to be accompanied by contemporaneous 
EUV/X-ray observations. 



5 



Ehrenreich et al.: Hint of a transiting extended atmosphere on 55 Cnc b 



55 Cancri e 



X 

_2 

"D 
CD 
03 

"CO 

E 
o 



o — 




U 2 



Chandrae-ray 



-4 -2 2 4 6 

Time from inferior conjunction (h) 

Fig. 8. Same as Fig. [5] for 55 Cnc e. Top: Lya light curve in- 
tegrated between 1 215.36 and 1 215.67 A. The vertical dashed 
line indicate mid-transit and the dotted lines the transit contacts. 
Bottom: Chandra count rate. 

4. Discussion 

4.1. Can 55 Cancri b not transit whilst hosting a (partially) 
transiting exosphere? 

We detect a Lya absorption signal during inferior conjunctions 
of 55 Cnc b, suggestive of a transit-like event. Can it be the sig- 
nature of a partial or grazing planetary transit with an impact 
parameter, b = a^lR* cos i, close to 1? If the orbital planes of 
55 Cnc planets are all alignecj^] and consid ering the incl i nation 
of 4 = 82.5+ \\ deg derived for 55 Cnc e by |Gillon etaL] ( |2012 i 
and the semi-major axis of 55 Cnc b (a/, = 0.1 1 au or lb. \\) «*), 
the latter could have an impact parameter in a range of [2.8, 4.0] 
(lcr) or [1.5, 5.2] (3(f), therefore excluding a transit for any rea- 
sonable values of the planetary radius. Meanwhile, owing to its 
Targe' separation with the star, 55 Cnc b also has a large Roche 
radius of ~ 1.2 Consequently, the Roche lobe of 55 Cnc b 
can undergo partial transits. If it is filled with hydrogen, as ob- 
served for hot Jupiters, then a transit signature may be observed 
at Lya. 

Besides, the value of i e found by Gillon et al.| ( |2012[ ) is the 
least favourable to a grazing transit of 55 Cnc b, assuming plan- 
ets b and e are aligned, among the other values found in the lit- 
terature ( |Winn et al.|201 1| |Demory et al.|2012) or in the recent 
report by the MOST team of an inclination i e = 85.4 ± 2.5 deg 
(D. Dragomir, 2012-08-27), allowing for an impact parameter of 
55 Cnc b lying between 1.0 and 3.2. 

We built a simple model to test this hypothesis. Considering 
the impact parameter b and the optically thick radius of the hy- 
drogen exosphere r exo as free parameters, we fit a transit light 
curve calculated with the routines of Mandel & Agol] ( p002| l to 
the data. A best-fit model with a^ 2 of 2.3 for 6 degrees of free- 
dom was obtained for b — 1.03 and r exo = 0.35 Roche radius. 
The result is shown in Fig. [9] 

Assuming 55 Cnc b has a radius similar to Jupiter, the pa- 
rameters given above would imply a small partial transit by the 



Mc Arthur et al. 



(2004) tentatively report an inclination of 53° ± 
for the farther out jovian planet 55 Cnc d from HST astrometry, 
indicating possible misalignments in the system. 
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Fig. 9. Lyman-a light curve of 55 Cnc b from Fig. |5](top panel) 
compared to a best-fit model transit light curve (thick line) with 
an impact parameter of 1.03 and an optically thick radius of 0.35 
times the size of the Roche lobe. Inset: Sketch of the system at 
inferior conjunction. The Roche radius of the planet is repre- 
sented by the dotted circle, while the shaded region corresponds 
to the optically thick hydrogen envelope. 



planet itself, that would be however challenging to confirmed or 
ruled out with high-precision photometry. Meanwhile, we em- 
phasise that this problem is underconstrained, since solutions 
with non-transiting planet b would still provide decent fits to the 
data. In addition, a few-degree misalignment between the orbital 
planes of these two planets cannot be ruled out. Confidently ex- 
cluding grazing transits of 55 Cnc b with precise photometric 
observations could bring a useful insight on this question. 

The simple model described above does not take into ac- 
count the dynamics of hydrogen atoms. Looking at the trans- 
mission spectrum from Fig. [6] the absorption signal takes place 
in the blue part of the Lya line, between and -76.5 km s _1 , as 
expected in the case of neutral hydrogen atoms repelled by the 
stellar radiation pressure towards the observer. In order to assess 
whether radiation pressure can reproduce the observed velocity 
distribution of the absorption signal, we calculated the 2-D tra- 
jectory of a hydrogen atom released from the upper atmosphere 
of 55 Cnc b with a tangential speed equal to the orbital veloc- 
ity of the planet. The particle is subject to the stellar gravity and 
radiation pressure. The latter is estimated by reconstructing the 
intrinsic stellar Lya line from the observed profile as shown in 
Fig. [6] following the method of |Ehrenreich et alT ( 201 1) . 

The ratio B between radiation pressure and gravity (Lagrange 
et al.||T998] > is reported as a function of radial velocity on the 
right axis of Fig.[6](top panel). Its value close to 1 between -80 
and +80 km s _1 and its steep decrease at larger velocities (com- 
bined with the planet's orbital velocity of 85 km s -1 ) implies 



that as soon as the atoms escape the planet, they are rapidly ac- 
celerated by radiation pressure at negative radial velocity down 
to -76.5 km s in ~ 120 h, as shown in Fig. 



10 



This accelera- 
tion time scale is just 2 to 3 times larger than the characteristic 
time scale of ionisation of neutral hydrogen by EUV stellar pho- 
tons, which is t ~ 50 h assuming a solar EUV flux, and we 
predict that escaping atoms must present a distribution of radial 
velocities ranging from to -76.5 km s _1 , in agreement with 



6 



Ehrenreich et al.: Hint of a transiting extended atmosphere on 55 Cnc b 




150 200 250 
Time (hours) 



400 



Fig. 10. Plot of the radial velocity of a hydrogen atom launched 
from 55 Cnc b in a simplified 2D simulation (solid blue line). 
The horizontal black dashed line shows the observed absorption 
limit velocity of ~ 76.5 km s _1 . The vertical blue dashed line 
shows the characteristic life time of a hydrogen atom for an ion- 
ising flux of 1 solar EUV flux. 



the limit of the detected absorption]^] Finally, we also note that 
the absence of any signal towards the red is consistent with the 
radiation pressure interpretation. 

We stress that additional data will be mandatory to con- 



firm this result. Since the star is variable in X-rays (Sect. 3.3 1, 
similar observations at a different epoch could result in an in- 
creased transit signal, or no signal, as seen for HD 189733b 
( Lecavelier des Etangs et aT1|2012| l, pleading for new simulta- 
neous Lya and X observations. If the transit signal of 55 Cnc b 
is confirmed, subsequent modelling efforts should help under- 
standing the depth and duration of the transit light curve as- 
suming atmospheric escape rate much lower than what is cur- 
rently known for hot Jupiters (see Sect. 4.2 1, taking into ac- 
count the roles of radiation pressure, ionisation, self-shielding, 
as well as stellar wind interactions as another possible acceler- 
ation mechanism of neutral H atoms at large distances from the 
planet (|Holmstrom et al. 2008 ; Ekenba ck et al.|20"T0 Tremblin 
& Chiang|2012) i~ 

In particular, the light curve tentatively hints at an early 
ingress, not unlike what was claimed to be seen during the UV 



transit of the hot Jupiter Wasp- 12b ( |Fossati et al.|20lG||Vidotto 
|et al.|2010| i. 



4.2. Implications of the non-detection of 55 Cnc e at Lyman a 

The absence of a Lya detection at the transit times of 55 Cnc e 
implies that if this planet has retained a thick atmosphere, its spa- 
tial extent has to be limited. Given its reported radius and mass, 



4.0+°J g cnr 



(Gillon 



55 Cnc e has a mean density of p, 

et al.|2012") . This is too low for this planet to be telluric, i.e., a 
scaled-up version of the Earth mainly consisting of an iron core 
and a rocky mantle. Instead, internal structure models presented 
in |Gillon et al. (2012 1 suggest that a large fraction (~ 20%) of 
volatile is needed to account for the mass and size of the planet. 
One possibility is that the planet is enshrouded in a thick enve- 
lope of supercritical (> 647 K) steam. This water envelope could 
then provide a large source of atomic hydrogen under the pho- 
todissociated effect of the incoming stellar EUV and X-ray ra- 
diation. In this scenario, inspired by the prediction of the obser- 



4 After a time /, there will be a fraction of exp(-?/r) of neutral atoms 
remaining. After 120 h, there will thus be ~ 10% of neutral atoms with 
a radial velocity of -76.5 km s -1 . 



vational signature of 'evolved oceans' on terrestrial exoplanets 
(Jura 2004| l, a Lya transit signature could be expected depending 
on the atmospheric mass loss rate and the stellar photodissociat- 
ing flux. 

We scaled the atmospheric escape model used in previous 



2003 



studies of hot Jupiters or Neptunes ( Vidal-Madjar et al 
|Ehrenreich et al.|201 1| |Lecavelier des Etan gs et al.|2012p to the 
orbital and physical properties of 55 Cnc e (a full description of 
the model will be provided in Bourrier et al., in preparation). Our 
non-detection allows us to set a 3-cr upper limit of S= 3 x 10 8 g s _1 
to the atmospheric mass loss rate of 55 Cnc e. 

Taking the stellar X-ray luminosity from Sect 



3.3 



we can 



calculate X-ray fluxes at 55 Cnc e and b respectively of 640 and 
12 erg s^'cm" 2 (in the 0.25-2.5 keV band). These fluxes cor- 
respond to energy-limited escape rates of ~ 5 X 10 s and 2 x 



10 7 g s -1 , respectively, using the prescription of Lecavelier des 
Etangs ( 2007 [ ) and |Ehrenreich & Desert| (2011 1. We note here 
that the energy-limited escape rate estimated for 55 Cnc b is 
much weaker than the rates of ~10 9 to 10 10 g s _1 observationally 
derived for the hot Jupiters HD 209458b or HD 1 89733b. 

For 55 Cnc e, our upper limit to mass loss is lower than the 
X-ray energy-limited rate, even without considering the EUV 



flux which may be an order of magnitude higher (e.g. Sanz- 
201 This appears to rule out high mass loss 



IForcada et al. 



efficiencies (<: 60%). 

A mass loss rate of 3 x 10 8 g s _1 is insufficient to completely 
erode a putative water envelope representing 20% of 55 Cnc e in 
mass. A rough yet conservative calculation, unphysically assum- 
ing that all the water (i.e., more than 1 M ffi ) has been photodis- 
sociated into hydrogen, yields a value of ^ 6 x 10 ll) g s for a 
complete atmospheric escape to occur within 5 Gyr. This is two 
order of magnitude larger than our 3-cr upper limit, indicating 
that a massive photodissociating water envelope scenario is not 
ruled out by our observations. At this stage, however, we note 
that a proper calculation would require a photochemical model 
for the upper atmosphere of a super-Earth such as 55 Cnc e. 

Finally, another interesting but more exotic scenario could be 
that of a thick carbon dioxide atmosphere on 55 Cnc e instead of 
a water envelope. In this case, we would not expect to observe a 
Lya signature during transits. 



5. Conclusion 

Our HST and Chandra observations reveal that 55 Cnc is bright 
at Lya and variable in X-rays. This opens interesting perspec- 
tives for testing the impact of the incident X-ray and UV radia- 
tion on the atmospheric heating and chemistry for planets in this 
system, which cover broad ranges of both mass and separation 
(see, e.g., the study of the GJ 876 system by Fra nce et al.|2012| l. 

No Lya transit signals are detected for 55 Cnc e, allowing us 
to estimate an upper limit to the escape rate of a gaseous enve- 
lope enshrouding this super-Earth. This limit, of ~3 x 10 8 g s _1 , 
does not preclude the existence of such a volatile envelope. 

We found that the Lya stellar emission is absorbed at 
times corresponding to the inferior conjunction of 55 Cnc b. 
Furthermore, the absorption is blueshifted and observed down 

to 76.5 km s _1 from the line centre, as expected from atoms 

expelled by the stellar radiation pressure. Since the planetary 
system is close to being seen edge-on, this hints at the tantaliz- 
ing possibility of a partially transiting extended H i atmosphere 
around 55 Cnc b, escaping the planet at a rate of 2 to 3 orders of 
magnitudes lower than what is known for hot Jupiters. 

Yet, additional observations at Lya, as well as in other lines 
(e.g., C n, Mg ii, or Na i) are clearly needed to confirm this result. 
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Meanwhile, it pushes to continue the photometric follow-up of 
55 Cnc, to firmly rule out grazing transits of 55 Cnc b in the 
visible, in case the orbital planes of planets b and e are slightly 
misaligned. 

If a transiting extended atmosphere on 55 Cnc b is confirmed, 
it would represent the first atmospheric characterisation of a gi- 
ant planet that is not a hot Jupiter: as such, it would contribute 
to bridge the gap between hot exoplanets and the cooler plan- 
ets in our Solar System, which experienced hydrodynamical at- 
mospheric escape early in their history. The bright star 55 Cnc 
would then become our best laboratory for these comparative 
studies. 
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